Functional Integration on Constrained Function 

Spaces 



J. LaChapelle 

1387 NW Ashley Dr., Albany, OR, 97321, USA 
Abstract 

Analogy with Bayesian inference is used to study constrained physical sys- 
tems within the context of functional integration. Since functional integrals 
probe function spaces, both kinematical and dynamical constraints are treated 
simultaneously and on equal footing. Following the analogy, functional coun- 
terparts of conditional and conjugate probability distributions are introduced 
for integrators and then applied to some well-known examples of constrained 
functional integrals. The analysis leads to some new functional integration tools 
and methods. These are utilized to construct a model of the prime counting 
function as a constrained gamma process. 



1 Introduction 



Kinematical constraints (that is, constraints in the form of boundary conditions) on 
physical systems modeled by differential equations have been well-studied. Indeed, 
from the degree of control that has been achieved, it hardly seems necessary to study 
them further. On the other hand, typical solution methods are based on elementary 
techniques that rely on simple boundary value matching or symmetry (as in the 
method of images). But subtleties can arise from more complicated geometries and/or 
topologies, so it is necessary to extend the elementary methods to include general 
geometries and topologies — especially in quantum physics. 

For example, one extension of elementary methods to general geometries makes 
use of the generalized Green's theorem: By formulating the solution of a differen- 
tial equation in terms of Green's functions, arbitrary boundary geometries with cer- 
tain regularity conditions can be treated. As another example, extension to general 
topologies that may introduce periodicity typically involves the study of functional 
determinants and spectral theory. In both cases these extensions deal directly with 
function spaces and the mathematical complexities and subtleties inherent in them. 

Similarly, dynamically constrained Hamiltonian systems and their quantization 
have been — and continue to be — extensively studied for obvious reasons. Solution 
methods for this constraint type are usually anything but elementary. The vast 
literature on this topic supports the contention that, here also, the function spaces 
of the dynamical variables (as opposed to their target manifolds) are of primary 
importance. 

From a function space perspective, the distinction between kinematical constraints 
and dynamical constraints is unnecessary. Both types can be formulated by posing 
a restricted or constrained function space: the restriction is often imposed indirectly 
on a target manifold, and it leads to some kind of set-reduction in some appropri- 
ate general function space. Consequently, one can anticipate that function spaces 
furnish a fruitful arena in which to formulate and study all constraints. Moreover, 
it has long been recognized that functional integration offers reliable if not always 
acceptably rigorous methods to study function spaces, so it is not surprising that 
functional integrals have become useful analytical and numerical tools to study com- 
plex, constrained physical systems. Importantly, they offer a means to incorporate 
both kinematical and dynamical constraints under one roof. 

There are many references in the physics literature that study constraints in func- 
tional integrals; largely utilizing heuristic approaches. The aim of this article is not 
to supplant those methods — they are useful tools — but to propose a mathematical 
basis for functional integration on constrained function spaces. The basis is suggested 
by analogy to Bayesian inference theory, and it affords some guiding principles. With 
guiding principles in place, useful integration methods can then be developed and 
tested against known results. 

This work will utilize the Cartier/DeWitt-Morette (CDM) scheme as the mathe- 
matical foundation of functional integration. A short summary is given in appendix 
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[A| but the reader is encouraged to consult p] [3] for background and details. Roughly 
stated, the CDM scheme uses Fourier duality to define linear integral operators on 
Banach spaces in terms of bona fide measures on Polish spaces. The as-defined func- 
tional integrals, which can be characterized by associated integrators, then inherit 
useful properties through their duality relationship. And these properties can be 
used to reliably manipulate the functional integralsjj 

Application of the CDM scheme to unbounded quantum mechanical (QM) systems 
is well-understood, but how it works under general constraints seems to require new 
principles. We start by presenting several well-known examples that contain clues 
to the underlying principles. To begin with, they suggest that constraints add non- 
dynamical degrees of freedom, and this requires an enlarged function space. Next, 
the Bayesian analogy suggests the notions of conditional functional integrals and 
conjugate integrators. Together with the functional counterpart of the Fubini theorem, 
these tools enable us to construct and manipulate functional integral representations 
of constrained QM systems within the CDM scheme. 

This is the main idea of the paper: Constrained QM systems require a state 
space comprising dynamical and non-dynamical degrees of freedom. Singling out a 
particular system by specifying particular constraints induces a subset of the general 
state space that we designate as a constrained function space. A specific QM system 
is then represented by a functional integral based on an appropriate integrator and 
constrained function space. Finally, the integral over the constrained function space 
is represented by a functional integral over the full state space but characterized by a 
conditional integrator by defining the functional analogs of marginal and conditional 
probability densities and using the Bayesian inference analogy to relate them. 

Needless to say, the proposed construction must reproduce known QM system 
results. So some implications and applications of the integral representations are 
checked against four archetypical classes of constrained QM systems. The four con- 
straint classes can roughly be characterized as i) kinematical, ii) dynamical, iii) peri- 
odic, and iv) discontinuous. In particular, we find it necessary to develop the notion 
of a gamma integrator which is a cousin of the Gaussian integrator. The proposed 
construction and its associated tools are then applied to particular QM systems. The 
subsequent efficient derivation of old results illustrates the utility of the new tech- 
niques introduced. 

As a new application of the formalism we conjecture functional integral represen- 
tations of some prime counting functions. The functional integrals represent average 
counting functions, and they give approximations to the exact counting functions 
that are superior to the conventional approximations. According to the construc- 
tion, the prime counting functions are modeled as a constrained gamma process - 

1 Since the CDM scheme is restricted to function spaces whose elements are pointed paths that 
take their values in some manifold X, i.e. maps x '. [t , t^] d M y X with a fixed point x(t a ) — x a £ X, 
the functional integrals in this article are strictly path integrals. However, the CDM scheme can be 
extended to include more general function spaces (see e.g. [4]) and the guiding principles we identify 
are not particular to path integrals in this restricted sense. So the term functional integral will be 
used interchangeably with path integral. 
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a perspective that may lead to a deeper understanding of the distribution of prime 
numbers. This is a nice realization of the often observed connection between physics 
and number theory (see e.g. [5]). 

Although the primary focus of the paper is a proposed construction of constrained 
functional integrals, there are several secondary results obtained along the way that 
we should point out: 1) The Gaussian and exponential integrators in the CDM scheme 
are redefined to include a boundary form and a parameter that encodes a mean path. 
The new definitions are more useful in the context of constraints, and they lead nat- 
urally to the notion of "effective action". 2) The complex counterpart of the new 
Gaussian integrator is defined. Although we do not pursue the idea here, it ap- 
pears that the complex Gaussian may contain important information regarding the 
Schrodinger^diffusion correspondence. Specifically, it might explain when analytic 
continuation succeeds or fails in this context. 3) In [6\ a new integrator within the 
CDM scheme was introduced based on analogy to a gamma probability distribution. 
Its utility for incorporating boundary conditions in path integrals was recognized, 
but its meaning and origin were unclear. Here we learn that the gamma integrator 
is a natural consequence of the Bayesian analogy. Moreover, the gamma integrator 
possesses a complex parameter that, when restricted to the natural real numbers, 
reduces to what can be characterized as a Poisson integrator. In consequence, the 
'propagator' for a QM system characterized by a gamma integrator yields an equiva- 
lent construction of the Poisson functional integral introduced in [I]. 4) Kinematical 
constraints suggest the notion of a Dirac delta functional on the topological dual of 
the constraint space. It turns out that these delta functionals are particular types of 
gamma functional integrals. In fact, gamma functional integrals can be used to define 
'step functionals' and all their (Gateaux) derivatives. Presumably, one could use such 
tools to construct a functional analog of distribution theory. 5) Dynamical constraints 
suggest the notion of Dirac integrators characterized by Dirac delta functionals on the 
constraint space (as opposed to the dual constraint space). Their heuristic equivalents 
have long been used to enforce functional constraints — the archetypical example is 
the Faddeev-Popov method. 6) QM systems with periodic boundary constraints can 
be efficiently encoded by paths in an appropriate fiber bundle that are constrained to 
be horizontal. 7) Discontinuous constraints can be handled using the idea of "path 
decomposition" ([7j-[9]). The functional integral tools lead to a recursive process to 
calculate propagators on segmented configuration spaces. The recursive process can 
be evaluated iteratively by the well-known Born approximation, but "path decom- 
position" suggests a new and potentially useful approximation technique rooted in 
boundary Green's functions. More importantly, the functional integral construction 
show that the propagators associated with such systems are intimately related to 
Poisson integrators and, hence, gamma integrators. 

(A caveat; all variables are assumed unit-less by appropriate normalization for 
convenience.) 
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2 Motivating examples 



2.1 Localization 

The first class of examples — constrained Feynman path integrals — can be charac- 
terized heuristically by the presence of a delta function(al) in the integrand of a path 
integral. Some particularly prevalent early examples in quantum mechanics of this 
type were point-to-point transition amplitudes, fixed energy transition amplitudes, 
and the propagator for a particle on S 1 . 

In the CDM scheme, the domain of integration for a path integral is a Banach 
space of pointed paths X a . So point-to-point transition amplitudes are obtained by a 
suitable delta function in the integrand that 'pins down' the loose end of the paths. 
Standard manipulations (pQ) show that the path integral for quadratic action can be 
expressed in terms of a restricted domain of paths 

/ 8{x{t b ),x b )e 2 ^ x '' x) Vu;{x) 

JX a 



where G(t&, t&) is the covariance associated with the gaussian integrator Vcu(x) de- 
fined on the space of point paths X a and W^ a ' b ^ is the variance associated with the 
space of point-to-point paths X ajb . 

Aside from the action phase factor and the resulting normalization!! 

[ Vj a ' b \x) = 1 , (2.2) 

Jx a , h y/det[iG{t b ,t b )\ 

the new gaussian integrator Vu ( - a ' b \x) defined on X ajb is characterized by a different 
covariance G {a ' b) {t,s) that exhibits the same boundary conditions as paths in X a ^. 

Now consider the other two examples. At the classical level, constraints such 
as fixed energy and paths on S 1 can be imposed by means of Lagrange multipliers 
in the classical action. It is then a standard heuristic argument that the Lagrange 
multiplier constitutes a non-dynamical, path-independent degree of freedom in the 
path integral that can therefore be integrated out. Essentially, this introduces what 
can be characterized as a Dirac delta functional. However, to give rigorous meaning 
to a delta functional, one would need a theory of distributions on X a . 

An alternative route is to define a Dirac integrator T>5(x) that does the duty of 
a delta functional. It can be thought of as a limit of a Gaussian integrator with 
vanishing variance, i.e. — > 0. The Dirac integrator reproduces the expected 
behavior; 

/ f(x)V5(x) = f(0) (2.3) 

JX a 

2 This particular normalization is fixed a priori from the choice f x T>ui(x) = 1. 



x„ 



-iriW (a ' b 1 (V) 



Jdet [iG (t b ,t b 



(2.1) 
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and 

where M : X a — > X a and M(x ) = 0, and it encodes a localization in the integral 
domain X a . 

Similarly, an inverse Dirac integrator T>5~ 1 (x) can be formally defined that corre- 
sponds to the case — > oo so that 

I e -^W^) VS -\ x ) = S(x') . (2.5) 

This integrator encodes a localization in the dual space X' a . In contrast to a Gaussian 
or Dirac integrator, this type of integrator is not translation invariant; 

I e -^'M)< D5 -^ x + x j = f e - 2 ™(*'+<>*)v5- 1 (x + x ) = 6{x') 

J X a ^ X a 

e- 2m{x °> x) V5-\x + x ) ~ V5-\x) . (2.6) 
Equivalent ly, 

[ e' 2m{x '^V5-\x + x ) = 5(x' - x' Q ) (2.7) 
Jx a 

where (x' + x' , x) = (x 1 , x + x Q ). 

The salient features of note from these three constrained Feynman path integral 
examples are: i) constraints are related to a localization in function space (or its dual), 
ii) constraints are related to a change in covariance and/or mean, and iii) in general 
the normalization of a constrained integrator is different than the unconstrained in- 
tegrator. 

2.2 Quotient spaces 

When the target space X of the pointed paths x : [t ,t&] — > X can be represented as 
the base space of a principal fiber bundle X, equivariant forms on X can be expressed 
in terms of associated forms on X. This technique is essentially symmetry based and 
allows determination of propagators on multiply connected spaces, orbifolds, compact 
Lie groups, and homogeneous spaces. 

The principal fiber bundle construction is essentially a generalization of the well- 
known method of images. In practice, the group structure allows the Banach space 
of^paths X a to be related to a group decomposition of the Banach space of paths 
X a where x : [t a ,tft] — > X. In this way, paths taking their values in configuration 
spaces with non-trivial topology and/or geometry can be treated as paths taking 
their values in the covering space. This renders a simplified function space — to the 
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extent allowed by the covering space. In terms of path integrals, the method can be 
roughly expressed as 

/ f(x)Vu(x) = I ff(x-g)Vu(x)Vg. (2.8) 

where G is the space of pointed paths g : [t a ,t&] — > G with g(t a ) = g a and G the 
group manifold. 

But the functions of interest are equivariant and covariantly constant which means 
X„ G X is parallel transported. Under this 'constraint', g t— > h G M.^ the holonomy 
group at x^, and the integral reduces to the standard result 

/ f(x)Vu(x) = [ [ f(x-h)Vu(x) dh . (2.9) 

JX a JW^ a JXa 

The point to be made here is that the two function spaces X a and X a are related 
through an integration (and/or summation for multiply connected or discrete holon- 
omy groups). Insofar as finite-dimensional integrals are localized functional integrals, 
we could loosely say that introducing a constrained integrator on X a x G along with 
some localization in G renders the constrained function space X a . 



2.3 Discontinuous spaces 

Our final class of examples is comprised of configuration spaces in which x(t) expe- 
riences some kind of discontinuity. Particular cases include bounded configuration 
spaces and barrier penetration. The previous two classes of examples gave only a 
vague hint of how constraints influence a path integral. However, this third class of 
examples yields valuable clues and insights. 

If we believe that a variational principle lies at the heart of the quantum— ^-classical 
reduction, then we should re-examine the variational problem in the context of con- 
straints. Consider a boundary in configuration space. For point-to-boundary paths, 
the correct formulation is a variational problem from a fixed initial point to a man- 
ifold in the dependent-independent variable space. This type of variational problem 
introduces a variable end-point that can be interpreted as a non-dynamical depen- 
dent variable that encodes the implicit constraints imposed by the configuration space 
discontinuities and boundaries. 

To formulate the variational principle for paths taking their values in a space X 
that intersects a boundary, consider the dim(n + 1) dependent-independent variable 
space N = XxR with a terminal manifold of dimension (n + 1) — k defined by some 
set of equations {Sfc(x, t) = 0} where k < n, x G X, and t G [t a ,t ]. Let 

I(x) = / F(t,x,x)dt 

Jt a 

be the functional to be analyzed. The extrema of I(x) solve the variational problem 
for point-to-boundary paths. In particular, for the case of X = M n , the variational 
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problem is solved by the usual Euler equations supplemented by 'transversality' con- 
ditions (see e.g. [TO]). 

There are two limiting cases of particular interest. When the terminal manifold 
in N coincides with a boundary (or surface) in X, then k = 1 and the transversality 
conditions reduce to 

F(t b , x(t b ),x(t b )) = -uVS{x{t b )) ■ x(t b ) (2.10) 

where v ^ is a constant and x(t b ) is on the boundary. For free motion, (I2.10p 
implies that critical paths intersect the boundary transversely. 

The other case of interest is when the manifold in the dependent-independent 
space is 'horizontal', i.e. x(t b ) is fixed and the terminal manifold is a line along the t 
direction. This clearly corresponds to a point-to-point transition between two fixed 
points contained in a bounded region. The terminal manifold is determined by k = n 
equations and the transversality conditions yield 

F{t b , x(t b ), x(t b )) = V 6 F(t 6 , x(t b ), x(t b )) ■ x(t b ) (2.11) 

where e is a unit vector in the x direction. If, in particular, F = L + E where L is 
the Lagrangian of an isolated physical system and E is a constant energy, then this 
is just the fixed energy constraint (dL/ dx l )x % — L = E. Consequently, the variational 
problem in this case is solved by paths with both end-points fixed that have fixed 
energy (pj). 

There are two lessons to learn from this: i) when boundaries are present, we will 
need to introduce a non-dynamical degree of freedom, and ii) the boundaries will alter 
certain expectation values of the paths according to the transversality conditions. 

At this point, the nature of the new degree of freedom is obscure. However, if one 
wants a functional integral to represent the solution of a second order partial differ- 
ential equation with non-trivial boundary conditions, then a consistent formulation 
emerges if one is willing to associate the new degree of freedom with a non-Gaussian 
integrator. It turns out that the new integrator is closely related to a gamma probabil- 
ity distribution in the same way that the Gaussian integrator is related to a Gaussian 
probability distribution. 

The nagging question is "Why a gamma integrator?" The examples have furnished 
some clues: not surprisingly they point to probability theory. If the answer can be 
understood, perhaps formulations of path integral representations of more general 
differential equations will become evident. 

3 Constraints as conditionals 

Consider a physical system with dynamical, topological, and/or geometrical con- 
straints. Postulate that such constraints introduce non-dynamical degrees of freedom. 
The obvious idea to incorporate these degrees of freedom in a functional integral con- 
text is to enlarge the functional space. Consequently, construct B = X x C a Banach 
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product space. The space X corresponds to what is usually thought of as the space 
of maps, and C will be a space of non-dynamical degrees of freedom induced by any 
constraints. In a probability context, this additional product structure would intro- 
duce conditional and marginal distributions. In our context, we expect analogous 
structures; about which little can be said until the nature of the integrators on C are 
understood. 



3.1 Probability analogy 

Here it is fruitful to develop an analogy with Bayesian inference theorylf) Momen- 
tarily pretend that B is a probability space. Let Qx(x) and 0c(c) be the marginal 
probability distributions on X and C respectively. Bayes' theorem implies 

FT i i \ ®c(c\x)O x (x) 
&x{x\c) = — 

f x ®c(c\x)<d x (x) dx 

=: C(c\x)Q x (x) (3.1) 

where Qx and 0c* are conditional probability distributions on B. This yields insight 
into the constraint induced normalization noticed in the examples. 

This result is not surprising, because a constraint could alternatively be formulated 
as a map M : X — > Y where y &Y automatically obeys the constraint. Then change 
of variable techniques can be used to show the two associated integrators are related 
through a functional determinant which is essentially C(c\x). This is standard, but 
it shows that the probability interpretation is consistent (at least with change of 
variables) and it lends credence to the analogy. 

So far, we have only made use of Bayes' theorem. To profit further from the 
analogy, consider an optical setup where plane monochromatic waves are focused onto 
an observation screen. We wish to study the nature of the light source by placing 
various non-conducting apertures between the source and the observation screen. 
Sooner or later we discover that under mild intensities the irradiance pattern on the 
observation screen is determined by the mean and covariance of the transmittance at 
each point in the aperture. Moreover, by changing the wavelength and/or intensity 
of the source, the resulting irradiance pattern can be predicted. 

The Bayesian inferential interpretation of these findings is that the conditional 
probability density or likelihood Qx(x\c) — which describes the irradiance pattern 
for a given aperture — can be factored as a product of a functional F(x) of the 
transmittance x and a conditional likelihood that only depends on the mean and 
covariance of the transmittance. In general, the statement is there exist sufficient 
statistics S s (x) such that 

§x(x\c) = F(x)Q^ ) (S s (x)\c) (3.2) 



3 The classic reference for the probability concepts introduced in this subsection is [IT]; especially 
Ch. 2. 
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where F(x) is a functional on X and Qs a (x){S s (x)\c) is a likelihood on S S (X) x C. An 
equivalent statement (by way of Bayes' theorem) is that the conditional probability 

density Qc(c\x) oc (c) @a^(x)(5's( a; )l c ) is a functional of sufficient statistics. 

There are two key pointgj illuminated by the analogy. The first point is the effects 
of a constraint can be inferred from a subset S S (X) C X given Qs s (X)(S s (x)\c) and the 
marginal probability density Gc(c). And the second is the marginal and conditional 
probability distributions on C belong to the same conjugate family, i.e. 

Qc(c\x) oc 9c (c) Q^)(S.(x)\c) . (3.3) 

There is great value in these two key points: Given a particular likelihood func- 
tion and a set of sufficient statistics, the conjugate distributions on C are highly 
constrained. In fact, consulting a table of conjugate priors for standard distributions, 
one can readily find the associated conjugate families. 

There are, no doubt, further lessons to be learned about constraints from the 
probability correspondence, but at this point we leave the analogy and return to the 
CDM scheme of functional integration (see appendix |A] for a summary). 

3.2 Conditional and Conjugate integrators 

Return to B a Banach space of pointed paths, and amend the CDM scheme with the 
definition^] 

Definition 3.1 Let B = X x Y be a Banach product space, and let each component 
Banach space be endowed with CDM scheme data. Define 

Q B (b,b>) Q B (b,b>) 

® X \ Y {x\y,x \y) := -— = 3.4 

Q Y[y,y) J x e B {b,b')Ve x ,z x x 

and 

Z B ,{b>) Z B ,{b>) 

Z X '\Y'[X \y):= = - — , , f-pr • (3.5) 

z Y'\y ) J x , Z B i{V) dnx<{x ) 

These two functionals define an associated conditional integrator by 

Qx\Y(x\y,x'\y')Ve x}YtZxlY x\y := Z X '\Y'(x'\y') . (3.6) 

B 

The space J- x \y{B) of conditional integrable functionals consists of functionals defined 
by 

F,(x\y):= [ & xlY (x\y,x'\y')d^x'\y')= [ ® B (6 ' d^x'\y') (3.7) 

Jb> JB>®Y(y,y') 



4 These points assume the system is not driven 'too hard' so that the probability distribution that 
characterizes the system does not change during trials or observations. 

5 The use of O in both the probability and functional integral context is meant to be suggestive, 
but it should be kept in mind that the same symbol is referring to two distinct objects that should 
not be confused. 



9 



where fi(x'\y') is a conditional measurqj on B' . And the linear integral operator on 
Fx\y{.B) is given by 

Fp(x\y)T>e x . Yt z x . Y x\y = / Z x ,\ Y ,(x'\y') d^{x'\y') . (3.8) 
b Jb 1 

Proposition 3.1 

J ®x\Y{.x\y,x'\y')V e ^ Y ,z xlY x\y = ^jyy J b Q s(b, b')V QB>ZB b 

Ve xlY ,z xlY x\y ~ 9y ^\^ 9 B ,z B & (3.9) 
In particular, since the integrator relation holds for any y' G Y' , 

J F n( x \y) v &x\Y,z X \Y x \y = J B F »( x \y) e Y{y,y')Ve B ,z B b (3.10) 

most often with (y',y) = S s (y) or (y',y) = for all y EY. 

Proof. The proof follows immediately from definition 13.11 and the relevant CDM 
definitions. □ 

Evidently expressing integrals like j B F(b)Vb when B is a product space requires 
knowledge of the 'marginal' and 'conditional' integrators on the component spaces. 
Of course, when elements in X and Y are independent, the conditional integrator on 
B reduces to a simple product of standard integrators on X and Y . But we anticipate 
that constraints induce a dependence between elements in X and Y . 

Now specialize to the case when Y represents non-dynamical degrees of freedom - 
perhaps due to constraints. As suggested by the optical diffraction example, postulate 
that the physical system is described by 'sufficient statistics^ and that 6y and ®y\x 
belong to the same conjugate family. Then knowledge of the 'likelihood' functional 
®x\Y(%\y, x'\y') implies knowledge of the conjugate family of Q Y (y, y'). Consequently, 
the heuristic integral J B F(b)T>b will be well defined in terms of known functionals. 

Accordingly, the Bayesian analogy suggests the definition: 

Definition 3.2 Conjugate integrators are characterized by 

®Y\x(y\x, •) oc <d Ss (x)\Y(S s (x)\y, •) <d Y (y, •) (3.11) 

where 

r e Y (y,y')V @Y ,z Y y = Z Y (y') (3.12) 



>Y 

and the proportionality is fixed by normalization. 

6 The conditional measure is well defined as the restriction of /i to the appropriate sub-cr-algebra 
over B' . 

7 In a functional integral context, 'sufficient statistics' is interpreted naturally as a localization in 
the space of paths precipitated by some constraint. 
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Note that this implies (by Bayes' theorem) 

9x\Y(x\y, •) oc B ((S s (y),x),-) . 
This property suggests the solution strategy 



(3.13) 




Fp(x)V Q - Z -x 



J B 




I F^y,x,-)e Y (y,-)1 ) e B ,z B b 



J B 




X 



x YJy 



G ll (S s (y),x)V eX:Zx x 



F^y, x , -)®Y(y, -)V 0YjZY y Vq X)Zx x 



(3.14) 



where the integral on the left is interpreted as a constrained functional integral, i.e. 
an integral over the constrained function space X, the third line employs functional 
Fubini (Prop. A. 3 in |6]), and G M can be interpreted as a constrained functional 
depending on the constrains only through sufficient statistics. This is the functional 
integral analog of (13.11) . Notice that the statement is equally valid with X Y if 
one knows some S S (X); hence suggesting an alternative solution strategy 



Both strategies will be employed in the next section. 

The important point worth emphasizing is that Qy\x and Gy belong to the same 
family of integrators, and they are simply related through the sufficient statistics 
that describe the integrator on X. This quickly narrows the search for an integrator 
associated with a particular constraint. 

However natural the probability analogy may seem, the usefulness of the defined 
functional integrals rests on their efficacy — which in turn depends on establishing 
physically relevant conditional integrators. We describe in detail four particularly 
prevalent integrator families in Appendix [B] and use them in the next section to 
re-examine the motivating examples of 

4 Examples revisited 

The value of the previous section is perhaps best appreciated by application to the 
familiar examples from QM outline in £j2j 




(3.15) 
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Eigenfunctions of position observables in quantum mechanics (for the semiclassical 
approximation) are characterized by their mean and covariance. When viewed as 
functions of time, they become paths in some configuration space. The semiclassical 
quantum mechanics can be described in terms of functional integrals with Gaussian 
integrators on the space of paths characterized by an assumed mean and covariance. 

Applying the previous section to the case Y a = C a , it is clear that the imposition 
of constraints in QM can change the quadratic form and/or the mean associated with 
an unconstrained Gaussian integrator. The notion of conjugate integrators enables 
one to construct the 'marginal' and 'conditional' integrators associated with those 
constraints, and hence make sense of constrained functional integrals. We will see 
that the type of conjugate family depends on the type of constraint. 



4.1 Localization 
4.1.1 Delta functional 

Definition 4.1 Let the space of L 2 ' 1 pointed paths x : [t a ,t&] — > X be endowed with 
a Gaussian integrator Vu>^q(x). A delta functional on C' a is a constraint that alters 
the covariance but not the mean of x G X a relative to T>uJx,q{x). (Here C' a is the 
topological dual Banach space of the space of pointed paths in the constraint function 
space C a .) 

The conjugate family associated with a normal distribution of known mean and un- 
known covariance is a gamma distribution, so the integrator associated with a delta 
functional is expected to be a gamma integrator j§ 

As motivation, consider the lower gamma integrator family. Put a — 1, and let 
L : T a R n with (/?', Lr) H> 2tuA • u and A = A*. Then formally, 

2W,oo(r) — ► / e- 2mX -"du = 5(X) . (4.1) 

with A G M. n and the integral over R™ understood as an inverse Fourier transform. 
On the other hand, 

■^"V=i£m = izm- (4 - 2) 

and the integrator V^ a ^i )00 (t) should be understood as a limit; 

/ T>Ja,p',oo(T) := Jim / Vja^^r) . (4.3) 
J |t s |-k» J 

So when tq is imaginary, T>^x,ip' ,oo{t) can be thought of as the functional analog of a 
two-sided Laplace transform implying 

£>7m/3',co(t) = lim . (4.4) 

Ta H \re\-+oo Det (if3') 



8 The gamma family is also conjugate for exponential-type distributions so the following analysis 
holds for action functionals that are not quadratic. 
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This justifies the definition: 

Definition 4.2 Let Det(zc') = 0. A delta functional on C' a is defined by 

S(c') := [ V ll>icl>00 (c) (4.5) 

and a Heaviside step functional by 



c a 



6 ( c ') ■= T^TTTT / ^7o,,c,oo(c) . (4.6) 

'Ca 



r(o) 



The delta functional vanishes unless (c', c) = for some c ^ 0, i.e. it picks out the 
zero modes of c' . The step functional vanishes unless |(c',c)| < c G for some fiducial 
c a G C. 

Remark that this definition suggests the characterization 



^(c'Ht^T / 2>7a,W,oo(c) 



(4.7) 



Evidently gamma integrators and their associated functional integrals could be used 
as a basis for a theory of what might be called 'distributionals'. 

Return to the (free) QM point-to-point propagator in R". Impose the constraint 
S(x(t b ) — x;,) where x& := x(t b ) with the choice 

- {P(x)',c) := 2m / [x(t b ) - x b ] ■ c(t) dt 

Jt a 

rib 

= 2iri [x(tb) — X{,] • / c(t) (it 

=: 2tu [x(t b ) - x b ] ■ c (4.8) 

where c(t) G M n and : X a — > C^. Clearly, the integral in definition 14.21 reduces to a 
simple delta function on K n in this case, and we can write 

Vco £tQM {x) := / e 27ri l a: ^- x '']^ c W*P 1Ao o(c)Pa; 2)Q (x) (4.9) 



where the right-hand side is to be interpreted as the integral of a conditional integrator 



on X a x C a . 



9 Since C' a is a polish space, the delta functional can be interpreted in a measure theoretic sense 
as the complex Borel measure associated with the identity element in the Banach space .F(C a ), i.e. 
Id = F M (c) = f c , e li > c '< c > dn{c') Vc e C a (see appendix). 
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'or s = i); 



JXaXCa 



[ [ ce 27Ti (( 5t ^ x - s) )'' c )vu SiQ (x)Vc 



J Ca J X a 



[ c [ e 27Ti (( St ^y' c )vco 0tQ (x)Vc 



J C a J Xq 



e «B(x(t 6 ))^ det [ lG (t bM )} f ce( mW{5t » Y ' c )Vc 



JC a 




(4.10) 



^det[zG(t 6 ,t 6 )] 



where we used functional Fubini, the boundary form in the fourth line appears because 
Q on the Hilbert space Xo with the new boundary conditions is not self-adjoint, the 
fifth line follows since Det(iW(5 tb )) 7^ 0, and x in the last line is the mean (actually 
critical in this case) path with boundary conditions x(t a ) = x a and x(t&) = x b . That 
is, 



It should be emphasized that G(tb,tb) is the covariance matrix associated with the 
original boundary conditions x(t a ) = x a and x(tb) = Also note that, since the 
mean path is non- vanishing, there is a phase associated with the boundary functional 
B (defined in Appendix [B]) that is not evident in the unbounded case represented by 



But this is still not the propagator. There may be more than one critical path, 
and in general there may be more than one mean path relative to an action functional 

10 Although this calculation has been done many times by time slicing, semi-classical, or linear 
mapping techniques; the point of repeating it here is to demonstrate that the as-defined functional 
integral tools allow the calculation to be carried out entirely at the level of the function space — the 
target manifold only makes an appearance through the boundary conditions imposed on the mean 
and covariance and the regularization/normalization of the functional determinant. 

11 Significantly, we did not have to expand about the critical path to do the calculation. Including 
the mean and boundary form in the definition of the integrator automatically handles this for us. 
But it does more. It tells us that, when Q — > S, the 'sufficient statistic' of import is not the critical 
path but the mean path. So, for example, the semi-classical expansion in terms of the mean path 
automatically accounts for self-interactions. In other words, once the mean is known, the Feynman 
diagram procedure (now without loop diagrams) is a way to estimate the moments of the integrator 
T>uix,s- As already remarked in Appendix [B] this is the essence of the effective action approach in 
quantum field theory. 





(t 6 - t a ) 
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S with the appropriate boundary conditions. Hence, for self-adjoint D on X a , 



e mB(x(t b )) 



K Q (x a ,x b ) = ======= . (4.12) 

The subscript Q has been included here to emphasize that the propagator is deter- 
mined by a covariance associated with specific boundary conditions, and it is a sum 
over x of all Gaussian integrators with x having the appropriate boundary conditions. 
This will be a recurring theme: propagators are represented by a sum over relevant 
parameters of an integrator family. 

The point-to-point propagator is a specialization of the more general integral 
(assuming convergence of the integral) 

K s (x a ,x b ):=Y, I S^{x)')Vu s>s {x) = I 8{P{x)')Y,Vu s , s {x) (4.13) 

2j J X a J X a g 

where Det /3(x)' = and the mean paths satisfy associated boundary conditions. It 
is to be understood as an integral over a conditional integrator on X a x C a . 

According to the functional Fubini, the order of integration can be interchanged 
for product integrators (the analog of independent joint distributions). Then, since 
the gamma family is conjugate for Gaussian likelihood functionals, the integral over 
X a will yield another gamma integrator — in which case the integral with respect to 
Vc is well defined. Explicitly, 

/ 6(P(x)')Vu Si q(x) = / Vu SjQ (x)V~f imx y j00 (c) 

J X a J B a 

=■ I (z~ mxY > c) )^Vli,oAc) (4.14) 



and the Gaussian expectation of exp{— i(/3(x)', c)} must lie in the family of gamma in- 
tegrators. Likewise, ( exp{— ir/sQ(x — x)}) must lie in the family of Gaussian 

integrators — which brings us to the next subsection. 



4.1.2 Uniform functional 

Definition 4.3 Let the space of L 2,1 pointed paths x : [t , t&] — > X be endowed with a 
Gaussian integrator Vcu Sj q(x) . A uniform functional on C' a is a constraint that alters 
the mean but not the covariance of x e X a relative to Vuj^q{x). 

The conjugate family associated with a normal distribution of known covariance and 
unknown mean is again a normal distribution so the integrator associated with a 
uniform functional constraint is expected to be a Gaussian. Essentially this means 
that the unknown mean of the conditional integrator on X a is normally distributed 
with respect to the marginal and conditional integrators on C a . 
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Definition 4.4 A Dirac integrator \ is defined by 

V5 5 (c) := Vu 5oo (c) =: 5(c - c)Vc such that lim Det[W 1/2 ]e^ w{c,) := 1 . 

(4.15) 



Proposition 4.1 The Dirac integrator is normalized j c V5 5 {c) = 1, translation 



in- 



variant T>5 5 (c — Co) = V5 5 (c), and furnishes the functional analog of a Dirac measure 

[ F^c)V5- c {c) = F^c) . (4.16) 

JCa 

Proof. The normalization is obvious and the translation invariance follows from 
the translation invariance of the primitive integrator. The third relation follows triv- 
ially from definitions; 



C' 



[ F„(c)V8- c (c) = [ Z(c')dti(c') 

1 dfi(c') 

= F,(c) (4.17) 
where the last line follows since c = c =>- e m ( c '^ c - c ) = \ \/c' G C' a . □ 

The justification for the term 'uniform functional' is obvious: Similar to the 
gamma integrator with a = 1 characterizing a delta functional on X' a , the Dirac 
integrator characterizes a uniform functional on X' a by J c VS s (c) = 1. 

More generally, the Dirac integrator yields 



(4-18) 

J ° a C ( c o) 



with M : C a — > C a a diffeomorphism and M(c — c) = 0. 

To see that Dirac integrators are conjugate integrators, use flB.201) in the context 
of sufficient statistics. To make this concrete we stipulate that 

= Dcx ■ S s (X a ) ->■ C' a ; G~l = D xc : C a -)■ S 8 (X a )' . (4.19) 

Insofar as the sufficient statistics characterize 'classical' observables in a QM context, 
this means the constraints are only correlated with classical observables. Now taking 
the G cc — » oo limit gives 

e -nQ c (c-m c}x ) _^ e -irQ(c-rn clx ) 5 ( c _ m ^ ^ ^ 



12 The Dirac integrator is improper in the sense that it is a limit of a Gaussian that requires 
regularization to achieve a sensible normalization. 
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where Q{c 1 ,c 2 ) = {{G cx D xx G xc )~ l c l ,c 2 ) =■ (G~ l ci,c 2 ) and = c + G cx D xx (x-x). 

Clearly, Q(c\,C2) = Qx{x\, x 2)\s B (x a ) because of f l4.19D so the exponential can be 
viewed likelihood functional 

Bx\c(S a (x)\c,-) = e-* Q *^ (4.21) 

where now x := {D xx G xc )~ l m c \ x . It follows that 

Qc\x(c\x,-)~ e -^ c - s) 5 5 (c) (4.22) 

and, hence, 

Vum^,Qc\x(c\x) « e-* Q ^-*)5- c {c)V{x,c) . (4.23) 

Since the mean is altered but not the covariance, the pair (Vu) m , t Q c . x (c),VcUc !00 (c)) 
are conjugate integrators and j c 5(c — c)Vc is a uniform functional constraint on X a . 
Clearly the Dirac integrator implements a dynamical constraint in the sense that it 
can be viewed as a Lagrange multiplier term in an action functional. 

An explicit example of this type of integral comes from the fixed-energy propagator 
on phase space often heuristically defined by 

G ps (q b ,q ai E) = [ S(H{q, P )-E)Vcj^ b \q,p) (4.24) 
J(Q,P) 

where T>co^ a ' b \q,p) is an appropriately defined point-to-point integrator on some phase 
space. 

Nothing is altered if we allow general action functionals. That is, we can define 
Dirac integrators as the limit S — > for exponential-type integrators. The quintessen- 
tial example of this type (which, however, is outside the scope of this article) is gauge 
fixing in quantum field theory. Assuming well-defined functional integrals for fields, 
a Dirac integrator is just the Faddeev-Popov 'trick' 

! V5-JM(g)) = V 1 —— (4.25) 



'Jo 



where M : G — > G and M(g — g) = 0. Of course, an admissible gauge fixing condition 
requires a single g Q , and it is standard to average the delta functional over g with 
respect to some (usually Gaussian) integrator. 

The whole Faddeev-Popov procedure can be readily interpreted from a conditional 
integrator viewpoint 

G(F(A))e iS{A) V(F(A)) = [ G(F(A))e iS{A) DetF'VA 

J A 

= [ G(c)e iS ^BetF' 5(c - F(A)) V(A, c) . 

JAxC 

(4.26) 
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4.2 Quotient spaces 

Let M be a connected manifold without boundary and G — > M — )• M a principal 
fiber bundle endowed with a connection. We wish to define the functional integral 
J~ Fji(rh)Vm where M a 3 m : [t a ,t&] — > M. The problem is that the base space is 
complicated in general: It may be very difficult or impossible to directly define an 
integrator on M. 

On the other hand, the covering space is usually easier to handle, and we assume 
that we can define an integrator so that the integral j M F^imfDm is well-defined. 
We also assume that F^(m)(i) furnishes a representation of G and is equivariant so 
that 

F At (m)(t)-g = p(g- 1 )F M (m)(t) (4.27) 

where p is a possibly non-faithful representation of G. 

Now, in the CDM scheme expressions like m(t) are shorthand for a parametrized 
curve: m(6)(t) = m a ■ E(t,6) where b G B a and S(t,6) : M — > M is a global trans- 
formation such that E(t , •) = m a . This parametrization allows integrals over the 
generally non-Banach space M a to be expressed as well-defined integrals over B a . 

The first point to make is that the parametrization is gleaned from the local 
structure of the bundle Uj x G where Uj C M and a local trivialization is given so that 
m(fo)(t) = (fh(b), g(b))(t). The parametrization for the first component is dictated 
by the manifold structure of M. The parametrization for the second component is 
fixed by requiring parallel transport of m(t) — since this will restrict paths to M if 
that is where they start. Consider an open set Uj C M, and let Ai denote the local 
gauge potential relative to the canonical local section Sj : M — > ML The equation for 
parallel transport, 

dft(t) = -s<(t)4(m(t))dt (4.28) 

clearly indicates the conditional relation between m and g, and it emphasizes the 
interplay between constraints and conditionals. 

To attack this problem we need to be more explicit about the parametrization of 
M a . 

Definition 4.5 Let {uji = 0} where u>j G AT*M a and i G {1, . . .,r} be an exterior 
differential system with integral manifold (B, M a ). This system defines a parametriza- 
tion P : B M a by 

P*Ui = Q\/i. (4.29) 

In particular, if B = X x Y and i — 2, the parametrization can be written locally 
on M as 

dm(x)(t) = X ( a) (m(x)(t))dx(t) a m(t a ) = m a , . 

dg(y)(t) = Y (b) (m(y)(t))dy(t) b g(t a , ) = g a { ^ U> 

where a G {1, . . . , m a }, b G {1, . . . , m b }, m a + m b < dimM, and the set {X( a ), Y^)} 
generates a vector sub-bundle V C TM of the tangent bundle. The solution of b4-3(J\ ) 
will be denoted m(x,y)(t) = m a • H(t,x,y) where E(t, x, y) : M — > M is a global 
transformation such that E(t a , ■, ■) = Id. 
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The parallel transport equation implies dg(t) ~ A-dm(t) which implicitly encodes 
the constraint through (I4.30p . This particular parametrization yields (within Uj) 

m(x)(t) = (m(x),g(m(x)))(t) = (m a ■ Z(t,x) , g a ■ Ve~ ™ dt ) (4.31) 

where E = 11(E). On the other hand, since m is a horizontal lift, 

m(a;)(t) = Si(m(x)(t)) ■ g(m(x))(t) . (4.32) 

It is clear that the restriction expressed by ( 14. 28ft implies g(fh(x))(t) G G V {t, x}, 

and m(x) : [t , t&] — > M(m a ) where H( ma ) — > M(m a ) A M is the holonomy bundle. 

We have learned that, given m(x) and a local trivialization, the constraint 5(M(g)) - 
5(g — g A - m) will lead to a path confined to an open neighborhood of a section that 
is isomorphic to the base space. And since the constraint only shifts the path along 
fibers, we expect it to be represented by a uniform functional which implies a Dirac 
integrator V5g e (M(g)) should be used. Moreover, M just effects translation on G a 
so the functional determinant of M' is trivial and the zero locus coincides with the 
holonomy group H( ma ). 

So, for a Gaussian integrator, 

Fji{m){t)Vu{m) ■= [ [ F^m(x))(t)5(M(g))VgVu^ Q (x) 

F^(si(m(x))(t) ■ KjVu^Qix) dh 



( m a) 



( m a) 



X a 

-1 



pih-^F^mix^Vu^Qix) dh (4.33) 



X a 



where m(x)(t) € M(m a ). 

In particular, the point-to-point propagator on M obtains for the familiar choice 
F M (m(x))(t b ) = 5(m(x)(t b ),m b ); 

K(m a ,m b ) := ld(m(t b ),m b )Vuj(m) 

J M a 

p(h _1 ) / 5(m(x)(t b ),m b )VoJx ! Q(x)dh 
Jx„ 



(ma) 



f P {Yi- l )K hb (m a ,m b )d\i (4.34) 



where the propagator K^ b (m a ,m b ) is the point-to-point propagator on M(m a ) - 
associated with the homotopy class of paths indexed by h^ — pulled back to M0 



13 It is legitimate to write Ky lb (m a , mj) instead of K(m a , mt) = K(si(m a ), Si(nT&) • hf,) because the 
horizontal lifting does not depend on the trivialization and, hence, the particular canonical section 
Si. So we are free to choose the trivial section. Recall that K(m a ,mb) represents the point-to-point 
propagator, so its evaluation may include a sum over relevant paths. 
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In most physical systems, the integrator T>Ux,q(x) is invariant under the restricted 
holonomy group H9 m s acting on X a . Then K(m a ,m.b ■ h°) = K(m a ,mf,) where h° G 
H° m -j and the propagator on the base space reduces to the well-known result (for 
X ■ G -+ C) 

^(rrr;,irr b ) = J2^ Ama (^m) = J2J2 dA ^ xAma ^ K ^^^ ( 4 - 35 ) 

Am a Am a ge& 

where A ma labels the representation of the monodromy group G = H( ma )/Kr9 ^ at 
the point m a and d\ ma its multiplicity. 

On the other hand, since the free point-to-point propagator on the group mani- 
fold is known to be exact, then by ( I4.12p we get a trace formula for point-to-point 
transitions 

E TTTV^n = EE^ B X A °-(^[h t ](5;,^) • (4-36) 

4.3 Discontinuous spaces 
4.3.1 Boundaries 

We wish to define an integrator for a space of pointed paths Mf with m : [t a , t&] — > M 
and <9M 7^ sufficiently regular. Experience from the previous sections indicates that 
we should take B a = X a x T a and consider a product manifold N = M x M. + . 

The task is to make sense of an integral of the form j Md F^{m)V>m. So, for a 
Gaussian integrator on X a , take a gamma integrator on T a to account for the boundary 
constraint. The parametrization of n := (m, r) : [t a ,t&] — )■ N can be written 



G?m(x)(t) = X( a )(m(x)(t))da;(t) a m(x)(t a ) = m a 
rfr(t) = F(r(t))dt r(t a ) = r a = ' 

and the integral defined by 

/ F^{m)Vm := / F M (n(x|T))8x|T(a:|T, -)Vq z x\t 
J ms Jb„, 

F^{n(x, r), •)Vu3^ Q {x)V'y a) p> iTB 

B a 



(4.37) 



/ z,Q 

=: /" Hp(n(S a {x),T))Vi at p tTg . (4.38) 

It remains to infer the nature of H ll (n{S s {x), r) and the associated integrator param- 
eters. 
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From the variational principle, we learned that the constraints impose transver- 
sality conditions on critical paths. We also learned that the physically meaningful 
path parameter is to be associated with r and not t. So the plan is to search for a 
conditional integrator for which0 

Qt\x(t\x, ■) cx e Ss(x)lT (S s (x)(t), •) e r (r, •) (4.39) 

where S S (X) is determined by mean paths. We replace critical with mean paths for 
the sufficient statistic because the quantum analog of the variational principle is 

with T defined in (SSM . 

To simplify matters, only the two limiting cases of transversal intersection and 
fixed energy discussed in $2] will be considered. Recall that these cases correspond to 
point-to-boundary and point-to-point paths respectively. Let n = (m, f) represent a 
mean path, and define the "first exit" time t 1 - by ^(t -1 ) := (hi^t- 1 ) where m(t ± ) = 
nie G <9M and t 1 - := T(t ± ). Recall the mean path m is also a critical path of Q(m(x)) 
since Q is quadratic, and there may be more than one critical path. 

The functional integral for a functional F® that takes its values on <9M is defined 
to be 

$ s (m a ):=/ F 9 ^ (m{x){r L )) Vu^^x) V 1gAoo {t) (4.41) 

JB a 

Likewise, for a functional Fp that takes its values in M\<9M, 

$\ 9 (m a ):= / F\ a {m{x){r))Vu m E, Q {x)V llfl . Td {T) (4.42) 

where Tq := t l . These definitions hold for Q — )■ S, but the mean paths are no 
longer critical. Of course, the suitability of these definitions rests on their ability to 
reproduce known results (see e.g. [B]). 

Like the quotient space example, the two in teg rals simplify considerably for prop- 
agators. For the point-to-boundary propagatoil 15 !. 

K a (m a ,m B ) = Af(m a ) / 5{m(x)(r ± ),m B )Vu 1fl ± tQ (x)V'y 0Aoo (T) 

JBa 

= A/"(m a )^/ -===== I>7o- w iB(*ft-L),oo(r) 

= ^K)E/ wvFn ^ d{lnr±) (4 - 43) 

Tf JR+ v / det2G(r ± ,r ± ) 



14 The same idea was implicit in the quotient space analysis. There, Qx\g was determined by 



lo 



The need for the normalization constant Af(m a ) can be established from dimensional analysis. 
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where the normalization constant J\f(m a ) = f g Kg(m. a , rn^dniB and the end-point 
is nie := m(x(r J -)). For example, if M C M", then the boundary term goes like 
B(m ) ~ \m.B — m a | 2 /r -1 . If there is more than one critical (or mean) path, care 
must be taken to split the integral over the boundary into regions associated with a 
particular critical path. 

For the point-to-point propagator, 

K(m a ,m b ) = / 5(m(x)(T b ),m b )Vu m B tQ (x)V'y 1) o tT g(T) 
JB a 

= / I A ,-n( =S V "fl-*iB{mV),T d {T) 

~ jT a v det * G{n,T b ) 

[•t 1 - TTiB(fh E ) 

= V / rfr 6 (4.44) 

^ B Jo v /deUG ! (rb,r fc ) 

where m^ := m(x(T b )). Incidently, for M C lR n , it is not hard to show using the 
explicit expression for B{fh E ) and Proposition 3.1 from [6] that the normal derivative 
of K evaluated on the boundary is just Kq as it should be. In essence, 

Vm a / T>^ 1 ^iB{m),T d {r) •n = A/'(m a ) / Z>7o,-™B(m),oo(T) • (4.45) 

JT a d JT a 

These point-to-point and point-to-boundary propagators are restricted in the sense 
that the paths are not allowed to penetrate the boundary. However, there are cases 
of interest when the boundary represents a discontinuity, and the paths are defined 
on both sides of the boundary. 



4.3.2 Segmented configuration space 

Often the target space of paths will have codimension-1 submanifolds that induce a 
decomposition of path space 0,.xi l) = X a in the sense that each Xa has its own 
integrator. Such is the case, for example, when a Gaussian integrator is defined in 
terms of an action functional and the potential in the action is discontinuous. More 
explicitly, x ( a ] G X { a i} is the pointed path x$ : [t a , t 6 ] ->■ (M«, m^) where M = [j . M« 
such that the intersection M (i) n = dM^ is a sub manifold of codimension-1. 

The objects of interest in this case are the propagators from the previous sub- 
section. Since it is known that the propagators are kernels for certain differential 
operators, Green's theorem provides a convenient starting point for the analysis. 

For an operator L defined in a bounded open region U C M acting on complex 
scalar functions from an appropriate function space, 

/ {L(j))Tp- [ <f>(L^p) = [ B((j>,lp) (4.46) 

where L* is the formal adjoint of L, and the functional form of B is determined from 
Stoke's theorem and the particular boundary conditions associated with the function 
space. 
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In particular, let Ui = UU^ be a bounded open region in R 3 with one surface 
§ = U 1 - 1 ) R of discontinuity. Choose ip to be the Green's function of V* in Ui with 
vanishing Dirichlet boundary conditions on <9Ui and (f> the Green's function of V in 
HJW with vanishing Dirichlet conditions on §. Then the theorem gives the Green's 
function Tp of V in U 1 - 1 ^ with non-vanishing boundary conditions on § 



jp = 0+ / -da (4.47) 
Js 

and in Ui \ U« 

jp= I ■ da . (4.48) 

This theorem has a simple interpretation in terms of functional integral repre- 
sentations of propagators: For bounded regions that allow paths to penetrate the 
boundary (e.g. the discontinuity we are dealing with) the total point-to-point prop- 
agator includes the restricted point-to-point propagator, which does not allow paths 
to leave the region, plus the potential on the bounding surface induced by all sources 
accessible to paths that are allowed to leave the region. This prescription is equivalent 
to the "path decomposition technique" used in [7], [8], [9]. 

In other words, Green's theorem can be used to partition the space of paths taking 
their values in a segmented configuration space into restricted and unrestricted sets. 
This is useful because the paths in the partitioned sets have particularly convenient 
boundary conditions and their associated propagators are relatively easy to calculate. 

Let be the restricted point-to-point propagators with Dirichlet boundary 

conditions on S, and the restricted point-to-boundary homogenous propagators 
for their respective regions UW. These are the propagators derived in the previous 
subsection. According to Green's theorem, the unrestricted point-to-point propagator 
from mi^ to m^ in Uj with one surface of discontinuity can be written 

<Vi J) ,m^) := ^^;(m«,m^) + f K^(m^\a)K^(a : m^) da (4.49) 

where my' G ILJW and a G S. Intuitively, the unrestricted point-to-point propagator 
within a bounded region IjW is implicitly determined by the restricted point-to-point 
and point-to-boundary propagators in U^. Similarly, the point-to-point propagator 
from to is implicitly determined by the restricted point-to-boundary propa- 
gator in U (i) . Note that Itfjp has non-trivial boundary conditions on §, but it still 
satisfies Dirichlet boundary conditions on <9U. 

Equation (14.49j) is a familiar expression, and it is often solved by way of a Born 
approximation. However, the path decomposition idea along with the fact that 
inside the integral is evaluated on the surface, suggests^ that we replace (inside 



16 



This argument is admittedly highly heuristic. 
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the integral) with where 




(4.50) 



where K^F) is the unrestricted propagator evaluated in U^, and 



(4.51) 



are evaluated in the region on the other side (->j) of m b . According to the path de- 
composition picture, r(a) and t(a) measure the (probability amplitude) contribution 
of pointed loops based at o that lie on either side of § and so|r| 2 + |t| 2 = l. 

Usually it is much simpler to find than to iterate (14.491) . For example, take 
the fixed energy propagator in R with a step potential V(x) = Vq6{x). Let k and ky 
be the wave vectors to the left and right of x = respectively. Then it is easy to see 
from the surface integral that t ~ 2^Jk$kv a / (&o + kv ) since this must be symmetric 
under k Q -B- kv ] and, hence, r ~ (k —kv ) / (ko+kv ) from the normalization condition. 

Now let U be divided into three regions. We can use the results for a single surface 
of discontinuity by covering U with two overlapping copies of Uj.. That is, each region 
contains only one surface of discontinuity. There are now six relevant propagators 
with appropriate boundary conditions. Their nature depends on whether or not the 
boundaries intersect <9U. Since we require Dirichlet boundary conditions for point- 
to-point transitions and if^ can only propagate across a single discontinuity, care 
must be taken to use the appropriate Ui for any given transition. 

To simplify, specialize to a planar geometry, and order the regions {1, 2, 3}. There 
will be two classes of propagators; half-space-type in regions 1 and 3 defined by 
Ui flUi', and unit-strip- type in region 2 defined by Ui n U v where Ui = U (1) U U (2) 
and U r = U (2) U U (3) . Choose the partition U = Ui U U w so that Ui = U U) contains 
the initial point x^) . By combining the single surface result appropriately, the point- 
to-point propagator for two surfaces of discontinuity can be written 



where is the restricted point-to-point propagator and K^J^ the restricted point- 
to-boundary propagator derived from (14.491) . 

It is important to remember that the propagator depends on all critical paths. 
For regions bounded by two planes, there is obviously a sum over all 'bounces' within 



K (D) (x {i) x (j 



S -K (D) (x (l) x 




(4.52) 
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the bounded region0. These bounce transitions are encoded in the K d ( j) {xa\a <y ^) 
propagator. Hence, poles of the convolution of relevant bounce propagators yield 
spectral information for their corresponding regions. 

In general then, the point-to-point propagator for U = US. is determined 
recursively from K^\ which is the standard unrestricted point-to-point elementary 
kernel in the region with vanishing Dirichlet boundary conditions on <9U, and 

K {D) (x {i) x U) ) — S- K (D) (x (i) x U) ) 

+ [ K^^)K^\x^)da^ 

(4.53) 

where the region U n _i is chosen to contain the initial point. 

There are cases when the iteration of (I4.49P can be summed explicitly. The method 
(see e.g. [12]) essentially boils down to finding a Poisson integrator that is valid 
everywhere in X a . To see this, let S(x) = Q(x) + V(x) describe some general (action) 
functional, and suppose the propagator Kq has been found in each HJW. Define 
V = Kq oVo Kq by its evaluation on ordered graphs in M, i.e. under the linear maps 
L n :T a —f R+ we have m(a;)(r) (m(x(r 1 )), . . . ,m(x(r„,))) =: (mi, . . . ,m n ) G M n . 
Then 

V{m(x(r))) K Q (m b , m k )V(m k )K Q (m k , m fc _i), . . . , V(mi)lifg(mi, m a ) (4.54) 

where m a < m 1; < ... < m^ are the time-ordered graph nodes. A time-slicing analysis 
([12]) when M = M. m shows the propagator has the same form given by 

K s (m a , m b ) := K s (r a , r h ) = (V) x , n (4.55) 

where (^a,^ is defined in subsection IB. 0.41 

But the time-slicing analysis is only straightforward in M. m ; otherwise there are 
well-known pitfalls. However, the target manifold independence of the right-hand 
side of (I4.55P suggests that the definition is correct for any manifold M. The point 
is, (I4.55P is more general than a perturbation expansion: it is defined in the function 
space rather than on the target manifold M. Consequently it offers potentially new 
insight and calculational techniques. 

In particular, an obvious conjecture is that Ks(r a , Tb) can be represented by 



K s (r a ,T b )~ / / y(r)2? 7a ^ i7l (r) da (4.56) 

for an appropriate contour CcC. This is a substantial generalization of the Poisson 
average because now the target space is C rather than R + . Not only do we get phase 



17 To the extent that the boundaries are exactly parallel and/or the planes extend to infinity, this 
gives an infinite sum which can be written analytically in the usual way as an inverse propagator. 
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information not carried by the Poisson integrators, but we don't have to restrict to 
C9',r> = A<M / ,r>. 

Here again is the thematic idea connecting propagators to summation over an inte- 
grator family. According to the conjecture, the analytic structure of Ks(m a , m&) could 
be encoded in the integral defined by Ksij a , r^; a) := T(l — a) f T V{r)V^ a ^^ Tb {T). 
Indeed, if V(t) = V G C, (/?', r) = A(Id',r), |Ar 6 | > 1, and 9£(Ar 6 ) > 0, we get 



K s {r a , n) 



r(i-a) 



c 



2ni 



[VV^ jTb (r)] da 



r(i oj) 

V I — - — : — 7(a, Arfe) da 



c 



V 



E 

_n=— oo 

oo 



2ni 

r(i + n)(-i)- 



-P{-n,\r b ) 



vJ2(-±) n P(nAn 



n=0 



V . 



(4.57) 



where the contour encloses the poles of j(cx, r&) on the negative real axis. And this 
agrees with definition (1B.48|) of the Poisson average (V)\ iTb . 

Parenthetically, it is interesting to note that the case |Ar&| < 1 can be handled in 
the same way by using the inverse gamma integrator, defined in the obvious way, to 
define an inverse Poisson propagator. Like the gamma integrator, the inverse gamma 
integrator is also a conjugate family for a Gaussian of known mean so everything 
hangs together. 



5 Prime examples 

Our goal in this section is to give functional integral representations of the prime 
counting function and some of its relatives. We formulate the counting functions 
in the spirit of quantum mechanical expectation values in the sense that they will 
represent the sum over constrained 'paths' with certain attributes. Specifically, the 
paths are conjectured to follow gamma rather than Gaussian statistics. 

5.1 Prime counting 

Postulate that the prime counting function is the expectation of a gamma process 
with unknown scaling parameter due to the constraint associated with counting only 
primes. The functional integral that enforces the constraint must be a gamma in- 
tegral because the conjugate prior of a gamma distribution with unknown scaling 
parameter is again a gamma distribution. Therefore, according to the general con- 
struction, the constrained functional integral that represents the expectation value 
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can be written as a constrained functional which is integrable with respect to two 
marginal gamma integrators. Similar to the QM point-to-point free propagator ex- 
ample, we put i((3(TbY,c) = i(rb — A(rt)) • c. Then let us define the average prime 
counting function to be 

fra / T ^7l,i/3(r 6 )',oo(c)r>7a-2-Id',r a ( T ) 

tr a / Z > 7a-1,-H' 1 7 1 0-) X>7l,i j 8(7i)',oo(c) (5.1) 

«/T a xC a 



7Ti(r 6 ) := tr Q (r) T6 := 



where (/3(t 6 )',c) = ((S tb ,(T d - A(t 9 )))',c) G R+ such that r a (t 6 ) = n and A(r 6 ) 
represents an unknown possibly non-homogenous scaling factor, and the a-trace is 
defined below. 

Using the definition of a delta functional, and recalling the treatment of K(r a , r&; a) 
from the previous section yields 

7Ti(r 6 ) = tr Q / (-1)°~ 1 7(« - 1, ~n) X > 7i,» J 9(7i)',oo(c) 
./c a 

= tr^-ir-^a-l.-Afo))] 

r(i-a) 



l) a_1 7(a - 1, -A(r 6 )) da 



/" 7rcsc[7r(a + 1)] 
= X +1 2^ + 1) ( -1 ) 7(q» -A(t;,)) rfa 

= -E i 7TT7-7(n,-A(r 6 )) 



^ (n)! 

n=l v ' 



= -E 



r(n) 

^T(n + 1) 



P(n,-A(r 6 )) 



(5.2) 



where the contour begins at oo above the real axis, circles the point {1} counter- 
clockwise, and returns to oo below the real axis; that is, it encloses the poles of 
r(l — a). The series converges absolutely since 



lim 

n— >co 



n\ | 7 (n + 1,-A(r 6 ))| 



(„ + !)! | 7 ( n ,-A(7i))| 





1 


= lim 


(ra + 1) 





Evidently, 



00 I f 

Tifo+i)-Tifo) = -£(^i)!y 

n=0 v / • / - 



A(r 6 ) = . (5.3) 



A(r 6 +1) i 

e-H n dt~—— - (5.4) 

-A(r 6 ) log(T 6 ) 



represents the density of primes at r&. But we cannot infer that its inverse is the 
expected gap between primes since this is a conditional statement. We will return to 
this issue later. 
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Notice the significance of choosing 0' = —Id' for the unconstrained gamma process: 
It forces r b — > —r b and it includes a phase of 7r with each event (which is associated 
with a pole). Roughly speaking, this calculation simply sums the positive integers 
appropriately constrained with a non-homogenous scaling factor and weighted by 
r(n)/r(n + 1) = l/n0 

A good and obvious initial choice for the scaling factor is A(r&) = log(r&). Interest- 
ingly, an even better approximation is given by the somewhat curious scaling factor 
A(t 6 ) := log(^) where 

oo 

n := J2(- 1 ) m ~ 1 ( mT ^ /m ■ ( 5 - 5 ) 

m=l 

Numerics show that the exact prime counting function 7r(r&) (apparently) randomly 
oscillates about 7Ti(r b ) within a bound of Itt-i (fh) — 7t(tj,) | < | li (-77,) — 7r(Tf,) | for almost 
all T b eR+ — at least up to r b ~ O(10 14 )H 

Unfortunately, we are unable to give an explanation (if there is one) of why log (17,) 
is better than log(r&). If such an explanation exists, it will likely be tied to pair 
correlations associated with the gamma process. 

Remark that it should be possible to express incomplete gamma functions as the 
eigenvalues of a projection operator on the space of L 2 (IR + ) functions with compact 
support with respect to a generalized Laguerre polynomial basis, and it might be 
profitable to think of n 1 (r b ) in that context. That is, 7t 1 (t^) would represent the prop- 
agator (as a sum over point-to-point transitions represented by incomplete gamma 
functions) associated with a non-trivial dynamical gamma process. Perhaps the origin 
of r b would be more evident from that perspective. 

Similarly, the average Chebyshev theta function can be represented by 



0{n) = tx a (\{r d )) n := tr a / A(r 6 ) £>7 a -2 -iOto) CO 



To 



Kn) I 2 ( ^ +2)) (-ir7K -A(r 6 )) da 



00 



E 

n=l 

00 

E 



C+2 

(n + l)! 

lQg(Tfc) 

J (n+l)! 



7 (n, -A(r 6 )) 

7(n, -log(r 6 )) . (5.6) 



18 Of course, these manipulations can be bypassed (although with some loss of motivation) by 
simply assuming the relevant Poisson average defined in terms of the Poisson integrator that will 
lead directly to the summation of incomplete gamma functions in the last line. 

19 The numerical calculations were performed using Mathematica 9.0 which doesn't support prime 
counting beyond this order. Nevertheless, the approximation 7Ti(fh) is rather impressive. 
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Consequently, 



2 m <T b 



2^ 2^ fw + lV 7(»,-log(r 6 / )) 
l/mlog(r 6 ) 



EE 

2 m <T b n=l 



fn + l)! 



■7(71, -1/m log(r 6 )) 



(5.7) 



Like 7Ti(r(,), these also provide excellent estimates for r& h-> 77,. Note that, since 
|log( r & )7( ra, -log(f ^))| < |log(r 6 ) 7(71, - log(r 6 ))| for all n e N, then 9{r b ) < 9(n) 
and ^(rfi) < ip(n). 

This approach can be applied to twin prime counting as well. Our reasoning is 
heuristic. We maintain the hypothesis that the occurrence of twin prime numbers is a 
constrained gamma process. But now, taking pairs whose difference is n = 2 will incur 
the twin prime constant C2 normalization according to the standard probabilistic 
argument. Additionally, we do not want to include the counting event between the 
primes so the trace should only include odd positive integers. So we should have a 
sum of the form 

00 

^nj~C 2 J2W 2 (n)j(2n-l,-\(T b )) (5.8) 

n=l 

where Wzfo) is an appropriate weight. From the prime counting case it is reasonable 
to guess that0 Win) ~ (—1)™/ (T(n + l)r(n + 3)), and a little numerical experimen- 
tation shows that the joint probability weight is well approximated by 

Win) = — y —i r . (5.9) 

v ; T(n + l)r(n + 3) v ; 

Continuing this heuristic for arbitrary prime doubles 7^(11,) within an interval 
2i < Tb — 2, the normalizing constant becomes 

p\i 

for prime numbers p > 2. Inside the sum, the key is to group all integers between 
the associated primes together. That is, we still sum over odd positive integers to 
exclude the collection of numbers between the two primes, and the weight factor for 

20 The phase factor (— 1)™ shows up because the phases of the two primes give a joint contribution 
of ( — l) 2n which no longer cancels the (—1)™ coming from the residue at the poles. 
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the joint probability stays the same. In general then, the prime double hypothesis is 



2(-l) n 



^I> + l)I> + 3) 
+ 2)] 



2cos[7r(^±i 



1 r(2±i + i)r(s±i + 3) 



7 (2n-l,-A(r 6 )) 
7(n,-A(r 6 )) , 



r b - 2 > 2z e N . 



(5.11) 



Note that only the normalizing constant depends on i, and the series converges ab- 
solutely for finite r b ; 



lim 

n— »oo 



T(n + l)r(ra + 3) 



T(n + 2)r(n + 4) 
However, since lim 



|7(2n,-A(r 6 ))| 



| 7 (2n 



-A(r 6 ))| 



lim 

n— >oo 



(n + l)(n + 3) 



A(r 6 ) 



rj) _ >00 7(2n — 1, — A(t&)) = — T(2n — 1), and the series 
2(_i)T(2n - 1) 



{a n } ■-- 



T(n + l)r(n + 3) 



= . 
(5.12) 

(5.13) 



does not converge to zero, then ^(u) diverges as r b — > oo. Moreover, although W(n) 
was essentiall y g uessed with some numerical justification, the exact weight W{n) (if 
it is different )o is certainly bounded by 



-iy 



T(n + 3)r(n + 3) 



< W(n) < 



T(n + l)T(n + l) 



(5.14) 



and the sums associated with both bounds diverge as r b — > oo. Therefore, given the 
hypothesis of the non-homogenous gamma process for locating joint prime numbers, 
we conclude that 

lim 7r 2 i(r b ) oo Vi G N (5.15) 

T(,— >-00 

for any W(n) between the bounds — including W(n). 

Owing to its probabilistic foundation, the prime double hypothesis cannot be 
confirmed unconditionally. However, given the success of the average prime counting 
function TTi(T b ), it appears likely that the hypothesis is correct. In particular, if we 
accept it for at least i = 1, then verification of the twin prime counting conjecture 
follows immediately — albeit conditionally — since the sum diverges with r b . 

Of course, one might argue that the hypothesis is just an alternative to the Hardy- 
Littlewood twin prime conjecture. However, the hypothesis is not asymptotic. With 
it we can statistically verify the Goldbach conjectural: 



21 The agreement between the expected and actual prime double counting numbers is good but 
not as impressive as for single primes. We conjecture that W(n) is indeed the correct weight and 
the degradation comes from fb due to pair-pair correlations. 

22 I acknowledge Jean-Pierre Zablit for bringing the Goldbach conjecture to my attention. 
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Theorem 5.1 If the occurrence of prime doubles is a non-homogenous gamma pro- 
cess, then every even number greater than 2 is asymptotically almost surely the sum 
of two primes. 

Sketch of Proof: Assume the contrary. Then there exists a 2r Q that is not the 
sum of two primes. Clearly, r D cannot be prime. Further, r G cannot be 'straddled' by 
a prime double (p, p + 2i) with p + i = t for some i e {1, . . . , r D — 1} since otherwise 
p+[p + 2i] = (t -i) + [(r G - i) + 2i) = 2t . 

But according to the constrained gamma process hypothesis, the probability den- 
sity of prime doubles straddling the point r Q is given by the absolutely converging 
series 

Prtn r ( ~ 1)W A + (n,r ) + A_(n,T ) 

P,(r ) = C 2i ^ r(n + 1)r(n + 3) ^ (5 - 16) 

n=l v ' v ' 

where 

A±(7i,r ):=[±7(2n-l,-A(r ±l))T7(2n-l,-A(r ))] . (5.17) 
So the expected number of prime doubles that straddle r a is given by 

To-l 



S(t ) :=X) 7r «(2r )P i (r ) . (5.18) 



8=1 



Now, for sufficiently large r D the expected number goes like S(t q ) ~ r /(log 4 (r )). 
Moreover, since r /(log 4 (r )) is monotonically increasing for sufficiently large r a , it 
only takes calculating S(t ) for a few small r G to see that S(r ) > 1 for all r > 6. 
Since the probability that r is straddled by at least one prime double is 1 — e~ s( - To \ 
we have a contradiction asymptotically almost surely. □ 

One can check explicitly up to some sufficiently large t& = 2r that the probability 
of a contradiction is essentially almost sure. For example at = 10 9 we find the 
expected number of straddling prime doubles S(rb) > 29000 where we used the under- 
estimate J2i C2i/i ~ 1 for simplicity. So the probability that the next even integer is 
not the sum of two primes is less than about 10~ 12500 . It is perhaps disconcerting that 
the conjecture cannot be settled with certainty by this argument, but it is comforting 
that the probability that it is false — beyond where one is willing to explicitly check 
- decreases exponentially like e~ CTfc ^ log ^' ri '^ 4 with c ~ 0(1) a positive constant 1^1 

One final implication: Since the probability associated with prime doubles only 
depends on the gap between them through C^j, the probability of twin primes at an 
interval [r b — 1, r b + l]/2 is the joint distribution to use for the conditional probability 
of two primes being separated by a gap. So the expected gap between prime p\ and 
Pi given p\ is P(pi + and it is easy to establish that P(p) _1 ~ log(p) 2 . Hence 
Cramer's conjecture is true on average. 



23 If the prime numbers are randomly distributed, then it is plausible that a proof of the Goldbach 
conjecture with certainty is not possible. This would help explain its resistance to attacks. 
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5.2 Bound estimates 

Note that for n > 0, 



1 1 
= 1 - -73= + E(- 1 ) m_1 (^) 1/m_1 > 1 " • ( 5 - 19 ) 

On the other hand, 

oo oo 

n = n + Y,(-l) m - 1 (mT b ) 1 / m < r 6 + v / 2^E(-l) m - 1 



m=2 m=2 

oo 



m=l 



= r b - ^ < r 6 + (5.20) 

where Cesaro summation was used. Hence, n ~ r b — 3y / 27^/4 and — r&| < v^t^ 
for r 6 > 2. 

For the Chebyshev functions with n > 2, 

0<— i = 1--— - 1 > 7 n,-log r b 

dr 6 log(r 6 ) V V fzf (n + 1)! 



1 / 1 

log(r 6 ) V 



< 1, (5.21) 

and 



2 m <r 6 



m V login 1 /™) V n 1/rn login 1 /™) 



= V — l - l- 

2 m <T b 

V log(r 6 ); 2 ^ m 

l/m 



1 - 



—)-y 



n 



\og(r b ) J n 2 £^ Tb m 



< I i - — — 111 ! 



< 1 . (5.22) 
For the prime counting function with n> 2, 



d-Kiin) n - 1 



dr 6 r b log(r 6 ) 



(5.23) 
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yielding 



Trifo) - tti(2) = li(r 6 ) - log (log (r 6 )) + log(log(2)) - li(2) . (5.24) 

< log (log ( n )) < log 1 / 2 ^). 



Consequently -Ki(r b ) < li(r 6 ) and 



li(r 6 ) - 7Ti(r 6> 



So for r h sufficiently large, we have the estimates 

0(n) - n < y/2r b 
i>{n) - n < \/2r b 



(5.25) 
(5.26) 



li(r 6 ) - m{n) < V^g(n) , \K(n) - *i(n)\ < V^g(n) (5.27) 



using the notation fin) := fin)- Loosening the bound to \fW b renders the estimate 
for ipin) valid for all r b > 2. 

These estimates enable comparison between the average and exact prime counting 
function: 

Proposition 5.1 Assume that prime counting is a non-homogenous gamma process. 
For n > 8, let 

(5.28) 

(5.29) 

(5.30) 
(5.31) 



and 

Then 
and 





M 






1 




71" 


log(r 6 ) 



7r(r 6 ) - Tii{T b ) < \n(r b ) - h(r 6 )| < a/ r 6 log(r 6 ) 



Win) - 7Ti(r 6 )| < 



asymptotically almost surely^ 



Proof: The variance for the trivial gamma process is equal to its mean. So the 
variance of the constrained proces^f] goes like 



a ~ Kiin) < li(r b ) < 



2r fe 



log(r 6 ) 

Then, for a confidence interval depending on N, 



^in) - niin) 



< No- < N* 



Vr 6 > 



2n 
log(r 6 ) 



(5.32) 



(5.33) 



24 Numerics appear to indicate even tighter bounds 
|7r(Tb) — 7t(ti,) | < My/n/ log(rb) for constant N and M. 



Tr(n) - Tn(n) < Ny/n/ log(n) and 

) — nin)\ < My/n/ login) for constant N and M. 
25 This is essentially the variance of the number of events for the corresponding Poisson process. 
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Hence 



|7t(t 6 ) - li(Tb) I < vr(r 6 ) - 7Ti(r 6 ) + 7Ti(t 6 ) - li(r & ) 
Similarly, 



< AT, 



2r h 



k(r 6 ) - 7r(r 6 )| < 7r(r 6 ) - 7Ti(r 6 ) 



+ 



log (7b) 



2r h 



'" l/ log(r b ) + log(r 6 ) 



(5.34) 



. (5.35) 



Since the upper bounds on N and M increase monotonically to oo, the confidence 
levels for the estimates increase asymptotically to 100%. □ 



This proposition has obvious implications for the distribution of primes in general 
and the Riemann hypothesis in particular. But much like the case of the Golbach 
conjecture, it only statistically confirms the hypothesis asymptotically almost surely. 
Again, if the prime numbers are randomly distributed, then perhaps one cannot hope 
to do better than this by relating zeroes of Riemann zeta to prime counting functions. 
However, Riemann zeta as a function need not be tied to prime counting, and it is 
reasonable to expect the Riemann hypothesis can be verified unconditionally and with 
certainty by other means. 



6 Conclusions 

Constrained physical systems were studied from a function space perspective using 
newly developed functional integration tools. The tools rely on the notions of condi- 
tional and conjugate integrators — the analogs of conditional and conjugate proba- 
bility distributions in Bayesian inference theory. These notions show the well-known 
Gaussian functional integrals to be only part of the picture: To describe constrained 
systems, one must be able to manipulate functional integrals over constrained function 
spaces with conjugate integrator families. 

Using the constrained functional integral concept, well-known QM results were 
re-derived efficiently at the functional level. Additionally, the framework allowed 
construction of a model for various counting functions associated with prime numbers 
that give improved numerical estimates and, hopefully, a basis for understanding 
prime distributions. The examples analyzed here point to the utility of gamma and 
Poisson integrator families, but it is likely that other probability distribution analogs 
will be useful. 

No attempt was made to develop methods to calculate non-trivial gamma func- 
tional integrals. That Z{r') is comprised of the incomplete gamma function and 
functional determinants and it is defined for complex parameters, points to broad 
complexity. Evidently the study of X>7 is an involved but important project. The 
perturbation expansion notwithstanding, the gamma functional integral can be ex- 
pected to yield new calculation techniques. 
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It would be fruitful to extend the concepts developed in this article beyond sim- 
ple QM. In particular, the domain of x and r can be altered in obvious ways to 
include quantum fields and loops. Similarly, the domain of X a can be extended to 
include matrix-valued functions — opening the door to matrix QM. Together with 
the complex Gaussian integrator and the complex nature of the gamma integrator, 
such extensions would appear to offer broad applicability and significant potential. 



A CDM scheme 

The CDM scheme (pQ, [2], [3]) defines functional integrals in terms of the data 
(B,Q,Z,T(B)). 

Here B is a separable (usually) infinite dimensional Banach space with a norm 
||6|| where h 6 B is an L 2 ' 1 map b : [t a ,tj G K. — > M with M an m-dimensional 
paracompact differentiable manifold. The dual Banach space B' 3 b' is a space of 
linear forms such that (&', 6)bGC with an induced norm given by 

H&'H = sup|(& , ,&>|/||&|| . 

Assume B' is separable. Then B' is Polish and consequently admits complex Borel 
measures fi. 

and Z are bounded, /z-integrable functionals : B x B' — > C and Z : B' — > C. 
The functional ©(&,-) can be thought of as the functional analog of a probability 
distribution function and Z(b') the associated characteristic functional. 

The final datum is the space of integrable functionals F{B) consisting of func- 
tionals F{b) defined relative to ji by 

F„(&):= J Q{b,b')dfi(b'). (A.l) 

B' 

If fi i— > is injective, then ^(B) is a Banach space endowed with a norm 
defined to be the total variation of /i. 

These data are used to define an integrator T>Q^b on B by 

j e(b,b')V e , z b:=Z(b'), (A.2) 

B 

This defines an integral operator J on the normed Banach space J-'iB); 

B 

j F^b) V e , z b := J Z(b') dfi(b') . (A.3) 

B B' 
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And the integral operator J is a bounded linear form on T{B) with 

B 



F,(b)V e , z b 



B 



< \\F„ 



(A.4) 



B Integrator families 

B.0.1 Gaussian family 

Definition B.l Let X a be the space of L 2 ' 1 pointed functions x : [t a , t&] C R — y X 
such that x(t a ) =: x a G X with X real. A Gaussian family of integrators Vux,q{x) is 
characterized M 26 l 



Q(x,x') = exp{2iri(x' , x)} 

Z s ,w{x r ) = ^et[W 1/2 }e 2ni{x '^ sW( ~ x,) (B.l) 

where (x',x) 6 I, s 6 ^ean patlE x has boundary conditions x(t a ) = x a 

and x(tfe) = x b , the functional determinant is assumed to be well-defined, and the 
variance 

W(x' v x' 2 ) = l - {(x[,Gx' 2 ) + (x' 2 ,Gx[)} =: (x',Gx') {h2} (B.2) 

where the covariance matrix G : X' a — > X a with 3?(sG) is non-negative definite. 
Associated with the variance is a precision 

Q(xi, x 2 ) = (Dx, x) {1;2 } - B(x 1: x 2 ) (B.3) 

with a boundary form 

B( Xl ,x 2 ) := (Dx u x 2 ) - (D^x 2 , x x ) (B.4) 

where D = G^ 1 and is its time-reversed adjoin^ relative to a Hilbert structure 
on X a . 

The integrator family is defined in terms of the primitive Gaussian integrator Dx ; 

Vu StQ {x) := e -Ws)Q(x-*)x>x (B.5) 



26 This definition uses a different normalization from the usual Gaussian integrator in the CDM 
scheme. Both definitions are valid: we choose this normalization because it seems more consistent 
with definitions of other integrator families and it highlights the role of the functional determinant. 

27 To remind; the mean path is also a critical path in this case since Q is quadratic. Also, most 
physics applications of functional integrals take place in unbounded configuration space where it 
is customary to put x = 0. However, this seemingly innocuous simplification cannot be taken for 
conditional functional integrals since x plays an essential role as we will see in 2) This could be 
anticipated from the analogy with sufficient statistics. 

28 The time parameter runs 'backwards' for because the order 1 — !► 2 has been reversed. For 
example, if D = d 2 /dt 2 , then B(xi,x 2 ) = 1/2 f i; 1 i 2 (ii = when x(t a ) = i(tf,) = 0. Otherwise 
B( Xl ,x 2 ) = l/2(x 1 x 2 \ t t b - ±ix 2 \ll) so B{x) = xx\\ b ^ 0. 
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where T>x is characterized by 

G(x,x') = exp{2ni(x',x) - (7r/s)Id(a;)} ; Z{x') = Vse-* 8 ^ . (B.6) 

The primitive integrator can be thought of as a Gaussian with zero mean and 
trivial covariance, i.e. VujQ t ia(x). Loosely, Vx is the integrator analog of the Lebesgue 
measure on R n . Note that W (and hence Det W), inherits the boundary conditions 
imposed on x, and note the normalizations 

f e -W*m*) Vx = ^ (B.7) 

and 

/ Vu s , Q (x) = v^Detfjy 1 / 2 ] . (B.8) 

JX a 

The symmetry between the functional form of Z(x') and the integrand motivates 
the standard practice in quantum field theory of defining the effective action functional 
r:I a xX^Cby 

- (vr/s)r(x, x') := 2iri(x' , x) - ttsW(x') . (B.9) 

Clearly the exponentiated effective action is nothing other than the characteristic 
functional of the Gaussian integrator. Notice that, if D is self-adjoint, T(x,x') = 
T(x,Dx) = Q(x,x). In particular, T(x) = Q{x) where x! := Dx, and the functional 
integral is easily evaluated once x is known. 

The desire to handle action functionals more general than quadratic functionals 
suggests defining more general exponential-type integrators in the usual way; 

Definition B.2 Let X a be the space of L 2,1 pointed functions x : [t a , t&] C K — > X 
such that x(t a ) =: x a G X with X real. An exponential family of integrators T>Ux,s{ x ) 
is characterized by 

Q(x,x') = exp{27ri(a/, x)} 

Z s , Ws (x') = ^et{Wl ,2 }e^ i{x '^ sWs{x,) (B.IO) 

where (x',x) G R, s G {1,0? ^ e mean path x has boundary conditions x(t a ) = x a 
and x(tb) = X6 ; and Ws{x') = (x', Gsx') where Gs is the connected Green's function. 
The integrator family is defined in terms of the primitive Gaussian integrator T>x; 

Vu s ,s(x) := e^ ls)s{x ~ x) Vx . (B.ll) 

The effective actior^ associated with this integrator family is 

- (7r/s)r 5 (x, x') := 2m{a/, x) - 7rsW 5 (a;') , (B.12) 

29 Remark that, evaluating at x' :— Dx for self-adjoint D, we have r,s(S) = S(x) . That is, the 
effective action has the same functional form as the action relative to T>Wx y s{x)- Consequently, the 
functional integral is readily evaluated if x is known. But x is hard to find for S. In quantum field 
theory, functional integrals are defined relative to T>u) Si q(x)] in which case x is easy to find but 
then r becomes a functional power series in S. In fact, the whole effective action technique is a 
realization of this trade-off — finding the moments of an unknown probability distribution associated 
with (xs, S) relative to a Gaussian distribution associated with (xq, Q). 
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and the inverse of Gs is defined by 



D s (t,t')Gs(t',s)dt l 



5 2 r*(x) 5 2 W s (x' 



Sx(t)5x(t') 5x'(t')5x'(s) 



— dt' = *(t,s). (B.13) 



To see how conditional Gaussian integrators work, form the product space X a x Y a . 
Suppose a Gaussian integrator on X a x Y a is characterized by a positive definite 
quadratic form Q with mean m and vanishing boundary term. Put m = (x, y) and 



G 



G X x C*xy 



TherS 

Q (0, y)-rh) = Q x (x - m x \ y ) + Q Y (y~ y) 
where Q Y (y 1^2) = (D yy y 1 ,y 2 ), 

Qx(xi,x 2 ) = ((G xx - G xy D yy G yx ) ; ./•,.. r 2 ) 

and 

m x \ y = x + G xy D yy (y - y) . 
So the Gaussian integrator on X a x Y a is 

Vu~${x,y) := e-^^-^V(x,y) . 

On the other hand, 

%,g y (!,):=e-W s )W») % . 
Therefore, the conditional Gaussian integrator is 

— e -{K/s)Qx(x-m xly )<p^ x ^y^ 



' DuJ rn x]y ,Q x \ Y {x\ 



which yields 



VuJ m ^ Qx ^(x\ 



-(n/s)Q a (m x]y )jy e j. 



(Qx + Qy) 



Q 



Y 



-1/2 



(B.14) 

(B.15) 

(B.16) 
(B.17) 

(B.18) 
(B.19) 

(B.20) 
(B.21) 



with Qa(m x \ y ) := T(m x \ y , Dm x \ y ). 

In particular, let M : X a — > Y a be a homeomorphism such that Qi = Q2 M. If 
Y a = X a then G xy = G yx = since the x are independent Gaussian variables. Also 
y = M(x). Then formally, 



(B.22) 



It can be shown that Qx and Qy are positive definite since G is positive definite. 
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Vu SjQl (x\y) e 



(n/s)Q 1 (x) 



3 (Vs)Q 2 (S) 



Det 



Q2 
Qi 



so we get the standard result for a change of covariance; 



-1/2' 



-(n/s)Q 2 (x-M(x)) ViX = 



02 

Qi 



-1/2 



(B.23) 



(B.24) 



where V>\X is the primitive integrator on X a . Obviously the same condition holds for 
1 -H- 2 with T> 2 x the primitive integrator on M(X a ). 

B.0.2 Complex Gaussian family 

The previous subsection took the parameter s G {1, i}@ This restriction can be lifted 
by defining a complex Gaussian integrator. 

Definition B.3 Let Z 2 be the space of L 2,1 pointed functions (z,z) : [t ,t&] C R — y 
M c such that (z,z)(t a ) =: (z a ,Zg) G M c with M c a complex manifold. A complex 
Gaussian family of integrators Vu^ qc^w) on W a = Z 2 is characterized by 

Q(w, w') = exp{2Tri(w' , w)} 

Z^ w c(w') = Det^^V™^"^™ (B.25) 

where w := (z,z) G W a , w' = (z',z!) G W' a , and (w',w) G C. The complexified 
variance W c (w / 1 ,w' 2 ) = (w[,G c w 2 ) where the complex covariance matrix G c has the 
block form 

r<c — ( C*zz Gzx 

with ^t((w[, G c w' 2 )) > and G c not necessarily Hermitian\ 33 \ As in the real case, put 



(B.26) 



Voo^ Q c(w) = e'^^-^Vw 
where T>w is characterized by 

®{w, w') = exp{2iri(w', w) - irld(w)} ; Z(w') = e^ 1 ^ 



(B.27) 



(B.28) 



At the level of functional integrals, evidently there is little difference between 
the real and complex Gaussian families. The value in the complex case comes when 
the domain of integration is localized yielding complex line integrals. The complex 
Gaussian can be extended to the more general complex exponential for S c exactly as 
in the real case. 



31 The ratio of the phase factors is obviously trivial since y — M(x). 

32 That Gaussian integrators based on non-negative definite real G can be defined for s € {1, i] 
reflects the validity of the Schrodingerodiffusion correspondence through analytic continuation. 
However, analytic continuation does not maintain this correspondence in general. It is natural to 
conjecture that the analytic continuation Schrodingerodiffusion correspondence will break down 
precisely when G zz and/or G ZZl defined below, do not vanish. 

33 If z = z* then (G c )t = G W . 
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B.0.3 Gamma family 

Definition B.4 Let T a be the space of L 2,1 pointed functions r : [t a ,t&] — > C C C 
such that r(t a ) = 0. Its dual space T' a is the space of linear forms £(T a ,C). A lower 
gamma family of integrators T>^ a ^^ d {r) is characterized byz3 

e(r,r') = exp{*(r',T>} ; Z a ^ T9 (r f ) = D J ( ^l^ )a (R29) 

where a G C, f3' G T' a , Tq := (/?', tq) for some fixed tq G T a , the functional determinant 
is assumed well-defined, and 7 (a, Tq) is the lower incomplete gamma function given 
by 

TM = ^ ^', + ir (B.30) 

n=0 v 7 

TTie integrator family is defined in terms of the primitive gamma integrator by 



V la ,p, T9 {T) :=r a e-^' T) Vr (B.31) 
where r a is defined point-wise by r a (t) = T(t) a and T>r is characterized by 

6(r, t') = e W {i(r', r) - (Id', r)} ; Z(r') = 1 . (B.32) 
An upper gamma family of integrators W a ^r )Tg (T) is defined similarly where 

r(a,7^) = T(a) -j(a,Td) (B.33) 
is the upper incomplete gamma function. 

The primitive gamma integrator Vt is just ^ , 7o,id',oo( r ) ( or equivalently £T w / )C ,(t)) ' 3r ' 
It is normalized up to a factor of T(0); 



r(o) J Ta 

but 



J t ^7o,^,oc(t) = 1 = ^ jT X>r ^,o(r) (B.34) 



J t ThfajrAr) = = J t OT^oCr) . (B.35) 

Proposition B.l Lei L G C(T a ,T a ) be a linear map and L* G £(T^,T^) zfo /or- 
mal adjoint. The primitive gamma integrator is invariant under r y Lt (loosely, 
V(Lt) = Vt). 



:',4 



This definition is somewhat modified from the original definition in [6]. 



35 Notice that f T £>7o.id', oo{ T ) = 1 is the functional analog of the formal, normalized integral 



r(5J Io° e " ^( mu ) : = !• in other words, Pt ~ cZu/T(0)u. 
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Proof. 



l r '^° f e l{T '> Lr) - {ld '> LT) V(LT) = [ e ^V,T>-<L*id',T>p (i/r 

JT a JT a 



To 

i(f',T>-(Id',T) 

:/:, 



Dt = 1 (B.36) 



Hence, 



□ 



/ F^Lt)V(Lt) = [ F^Lt)Vt . (B.37) 

JT„. JT a 



Evidently, whereas the T>x integrator is the infinite dimensional analog of the 
translation invariant measure on R™, the Dt integrator is the analog of the scale 
invariant measure on Wl defined to be the identity component of GL(n, M). This 
scale invariance is responsible for the prominent role played by the gamma integra- 
tor in solving partial differential equations. Or, reversing the interpretation, the 
gamma integrator (which is dictated by conjugacy) is responsible for the asymptotic 
'time' behavior of parabolic and (bounded) elliptic/hyperbolic second order differen- 
tial equations. 

Put B a = X a xT a . For 0^ Gaussian and 6^ gamma, use the relation for conjugate 
integrators to get 

/ & X \t(x\t, -)Ve x]T ,z x]T x\T = / e x|T (x|r,-) — — —V eB ^ B b 

JB a JB a Z T' 

oc / Qs 3 (T)\x{Ss{r)\x,-) V @BtZB b. 

JB a Z T' 

(B.38) 

This suggests that integrals of conditional functionals on X a x T a be understood as 
f F^x\t)V @XIt>Zx]t x\t = / F^(S s (r), x)Vco^ Q (x)V la ^^(r) (B.39) 

J B a J B a 

when S s (t) is a sufficient statistic for the integrator family characterized by Qx- This 
is just a specialization of the solution strategy (I3.14p . and it plays a prominent role 
in the solution of differential equations. 

B.0.4 Poisson family 

Take the lower gamma integrator and regularize by replacing j(a, tq) with the regu- 
larized lower incomplete gamma function P(a,Tg) := 7(0, tq)/T(oc). Restrict to the 
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case a = n G N, f3' = Id', and tq = u G K+. Note that, for iV G Poisfa) a Poisson 
random variable, we have 

Pr(iV < n) = e ~ n ^kT ■ ( B - 40 ) 

k<n 

Hence, 

f 1 ^ 1 /* 

Pr(iV > n) = J2 e kT = P ^ = f(n) ^ 

k=n ' ' \ > Jt » 

which, in particular, implies 

a fc=0 

On the other hand, 

^-w = -*rL -L dTl '-' dTt - (B43) 

Not surprisingly, Pois^) is closely related to the restricted gamma integrator 
which motivates the following definition: 

Definition B.5 Let T a be the space of L 2,1 pointed functions r : [t a , t^] — > R + := 
M + U {0} endowed with a gamma integrator with a = n G N and (/3',r) = A(Id',r) 
such that A G C + @ The Poisson integrator V7i n x, Tb ( T ) is characterized by 

0(r,r') = exp{*(r',T>} ; Z n , x , Tb {r>) = P (n / ^ . (B.44) 

Det (Id — ^r'l 

TTie Poisson integrator is defined in terms of the primitive gamma integrator by 

-Dn nXTb (r) := We^'^Vr . (B.45) 

Note the normalization 

f Vn 0XTb (r) = l (B.46) 

and 



f V7r nXTb (r)=P(n,Xr b ) . (B.47) 

Given a functional F^(r) integrable with respect to £>r, define its Poisson average 
(i^by 

oo „ 

(F) XiTb := ^(r) P7r n , Vi ,(r) . (B.48) 

n=0 



:i(> 



C+ := M+ x iM is the right-half complex plane, i.e. zeC such that Re(z) > 0. 
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Now let L n : T a -> R™ by r m- t := {n,...,r n } where r< = (r/,r) = (L*(Id'),r) 
such that < T\ <,...,< r n < r fe and ^™ =1 7f = r&. Then under L n , F^{t) \-t 
i^(L n (r)) =: F(t), and using obvious notation for the measure in R™ 

r I \\ n \ n r 

/ F M (L n (r))P7r„, AiTi) (r) -> e~ A ^ 1 \ _6(r b - t)F(t) dr (B.49) 

where the symmetry factor (— l) n /n! comes from the phase e 7 ™ (associated with the 
underlying complex nature of r(t)) and the counting factor associated with inter- 
changing components of r. 

So, if F M (r) happens to be such that F M (L n (r)) is given for all n, then 

(F) x , n -> £ ^ [_B(r b - r)F(r) rfr . (B.50) 

ill I mn 

n=0 ^ M + 

In particular, if F m (t) = 1, 

°° \n /■ °° n \n 

(J^-ye-^V;^ /_e(T 6 -T)dr = e-^y;^- = P(0,A7i) = l. (B.51) 

n=0 ,/JK + n=0 

On the other hand, 

oo „ oo 

(Jd>A, n = E^- 1 )" / ^A,. b (r) = ^(-l)»P(n,A7i) . (B.52) 

n=0 n=0 

There is no inconsistency here because 

oo 

^(-l) n P(n,Ar 6 ) = [P(0,Ar 6 )-P(l,Ar 6 )] + [P(2,Ar 6 )-P(3,Ar 6 )] + --- 

(Ar 6 ) 



n=0 

00 " \2n 



E 



2n! 

2n=0 

= P(0,Ar 6 ). (B.53) 

This definition of a Poisson functional integral agrees with the definition in [1] 
and, hence, gives an alternative characterization of a Poisson integrator. It is already 
known that such functional integrals give solutions to first-order operator differential 
equations. The fact that the Poisson integrator is composed within the family of 
gamma integrators indicates that gamma integrators will play a role in solving first- 
order operator differential equations with constraints. 
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